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ABSTRACT 

We study the variability of major atmospheric absorption features in the disk-integrated spectra of 
the Earth with future apphcation to Earth-analogs in mind, concentrating on the diurnal timescale. 
We first analyze observations of the Earth provided by the EPOXI mission, and find 5-20% fractional 
variation of the absorption depths of H2O and O2 bands, two molecules that have major signatures in 
the observed range. From a correlation analysis with the cloud map data from the Earth Observing 
Satellite (EOS), we find that their variation pattern is primarily due to the uneven cloud cover 
distribution. In order to account for the observed variation quantitatively, we consider a simple opaque 
cloud model, which assumes that the clouds totally block the spectral influence of the atmosphere 
below the cloud layer, equivalent to assuming that the incident light is completely scattered at the 
cloud top level. The model is reasonably successful, and reproduces the EPOXI data from the pixel- 
level EOS cloud/ water vapor data. A difference in the diurnal variability patterns of H2O and O2 
bands is ascribed to the differing vertical and horizontal distribution of those molecular species in the 
atmosphere. On the Earth, the inhomogeneous distribution of atmospheric water vapor is due to the 
existence of its exchange with liquid and solid phases of H2O on the planet's surface on a timescale 
short compared to atmospheric mixing times. If such differences in variability patterns were detected 
in spectra of Earth-analogs, it would provide the information on the inhomogeneous composition of 
their atmospheres. 

Subject headings: Earth - scattering - techniques: spectroscopic 



1. INTRODUCTION 

Determining the nature of the atmospheres and sur- 
faces of exoplanets is of primary importance in prob- 
ing not only their formation history but also in identi- 
fying possible signatures of life. While it is very chal- 
lenging, and perhaps only feasible through direct photo- 
metric/spectroscopic observations of the planetary light 
resolved from that of the host star, there are several pro- 
posals for eventual astrobiological investig ations of po- 
tentially habitabl e, rocky exoplanets (e.g., Levine et al.l 
[2009: Savranskv et al.ll2QlQ : Matsuo fc Tam ura 2010l). 

The available exoplanetary light would be disk- 
integrated, i.e., that of a point source without any spatial 
resolution. Deciphering the exoplanetary light properly, 
therefore, is inevitably a highly difficult task, in partic- 
ular for those planets with diverse surface types and at- 
mospheres like our own Earth. 

Indeed, habitable planets are likely to exhibit a variety 
of complex patterns of their surfaces and atmospheres, 
intrinsically dependent on their climatology and geol- 
ogy. According to the global water cycle, liquid water 
on the planetary surface vaporizes, forms clouds, is car- 
ried by atmospheric circulation, and precipitates as rain- 
fall/snowfall. Depending on the total amount of water 
and the atmospheric circulation pattern, the surface of 
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the habitable planet s may be partially covered by ocean 
(e.g. A be et al.|[2QT l). and be observed only through at- 
mospheres with highly inhomogeneous and variable cloud 
cover patterns. 

Given these complexities, techniques to properly de- 
cipher the disk-integrated light of exoplanets need 
to be developed. Towards that goal, the time 
variation of planetary light due to spin rotation 
and orbital revolution is a powerful tool, and sev- 
eral authors have computed the expected variation 
patterns and pr o posed the r econstruction methods 
(iFord et al.l 120011: i Tinetti et al.' '2006a. b; Cowan et al] 
20091: lOaklev fc Cash 2009; Kawahara Fuiii 2oTo[ 
Robinson et al.ll2010l: Cowan et al.ll201ll: iRobinson et al.l 



20111 ; iKawahara fc Fuiii j 120111: iFuiii et all l2010L ijuiT 
Fuiii fc Kawahara 2012: ISanroma fc Pallell2012D . 

The variation of the continuum level in the visible to 
near-infrared (NIK) range mainly reflects the distribu- 
tion of landmass, ocean, cloud cover, and possibly vege- 
tation. The peak-to-trough diurnal variability for 0.1/im- 
wide photometry of the Earth in the visible / NIR range 
is found to be 10-30% (Livengood et al.""201l). Thermal 
emission of the Earth also shows a few percent of diurna l 
variation in the mid-infrared ( Gomez-Leal et ani2012f ). 
which primarily originates from the uneven cloud cover 
and humidity. 

In addition to the light-curve in broad-band photome- 
try mentioned above, molecular absorption depths ex- 
hibit diurnal variation. Since atmospheric absorption 
depths of molecules are determined by the column den- 
sity of the corresponding molecules along the optical 
path, they are sensitive to the presence of highly reflec- 
tive cloud cover in the atmosphere that effectively blocks 
the spectral influence of molecules in the lower atmo- 
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sphere. For idealized uniformly mixed atmospheres, the 
diurnal variation of molecular absorption features should 
be strongly linked to the spatial distribution of clouds. 

In reality, the distribution of molecules itself is not 
entirely uniform and may be time-dependent. Moreover, 
the spatial and time variations are not the same for differ- 
ent molecular species comprising the atmosphere. In par- 
ticular, the local column density of water vapor is known 
to va ry widely on a timescale of ~ 1 hr (Blake & Sha^ 
l2Qllf ). while other molecules, such as N2 and O2, are 
well-mixed in the troposphere and thus show no signifi- 
cant variations on short timescales. 

This paper extends our previous wo rk on photomet- 
ric lisht-curves in visible/NIR bands (iFuiii et al.ll2QlQL 
[Mil ), and examines the diurnal variation of molecu- 
lar absorption signatures in the disk-integrated reflection 
spectra of the Earth with future application to Earth- 
analogs in mind^. We focus on the absorption bands 
of H2O at 1.13/im and O2 at 1.27/im, which are among 
the most prominent absorption bands in NIR and are an 
indicator of habit ability and a biosignature molecule, re- 
spectively. We first analyze data from space-based NIR 
spectroscopy of the Earth by NASA's EPOXI mission. 
Then we consider a simple model (referred to as opaque 
cloud model) that reasonably reproduces the observed 
variation pattern, if the cloud pattern data are provided 
separately. Our model implies that the different behav- 
ior of water vapor and oxygen in their absorption band 
variations can be ascribed to their intrinsically different 
spatial distribution patterns, rather than to the common 
cloud coverage. We discuss how this difference can be 
used to probe signatures of surface and atmospheric in- 
homogeneity of exoplanets with next-generation direct 
imaging. 

The organization of this paper is as follows. Section [2] 
introduces the EPOXI data used in this paper and analy- 
ses the correlation between absorption depths and pixel- 
to-pixel climatological data obtained with Earth Observ- 
ing Satellites. Section [3] describes our model to reproduce 
the variation pattern of absorption depths and compares 
the simulation results with observation. Section |4] fur- 
ther discusses the different behavior between H2O varia- 
tion and O2 variation. Finally, Section [5] draws our main 
conclusions and discusses the implication for future ob- 
servation of exoplanets. An analysis of CO2, a molecule 
with properties intermediate between those of H2O and 
O2 in some ways, is presented in Appendix [Al 

2. ANALYSIS OF EPOXI DATA 

2.1. EPOXI observation 

The EPOXI^ mission (Livengood et al." '2011") per- 
formed photometric and spectroscopic monitoring of the 
Earth and the Moon from space as a benchmark for fu- 
ture characterization of terrestrial exoplanets. A part 
of the observations were devoted to spectroscopy of the 

^ The term "Earth-like" has been used extensively in the litera- 
ture, often without any precise definition. In this paper, the term 
"Earth-analogs" is used to refer to rocky planets that resemble 
the Earth in having surface temperatures, obliquities, continents, 
oceans and atmospheres sufficiently like the Earth's to give them 
similar global hydrologies. 

^ The Deep Impact flyby spacecraft for the Extrasolar Planetary 
Observation and Characterization investigation (EPOCh) and the 
Deep Impact extended Investigation (DIXI) 
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Fig. 1. — Examples of NIR reflection spectra of the Earth 
observed by EPOXI in March of 2008 (Earth 1: equinox) 
(Liveng ood et aD l2011': 'Rob inson et al .11201 11 ) . Different lines rep- 
resent different short exposures obtained at 2-hour intervals and 
labeled by the time tlhr] from the start of the observations on 
that date (see Figure [3]). Three vertical dashed lines indicate the 
wavelength ranges that we adopt to compute the equivalent widths 
(H2O: 1.07-1. 24/im, O2: 1.24/im-1.283/im, CO2: 1.59-1.62/im). 

disk-integrated scattered light of the Earth over the 
wavelength of 1.10-4.54/im. These observations were car- 
ried out on 2008 March 18-19, 2008 June 4-5, 2009 March 
27-28, and 2009 October 4-5, with 12 exposures every two 
hours of the day (integration time for each exposure is 
less than 2 sec). The sub-observer latitudes (the latitude 
of the intersection between the Earth's surface and the 
line connecting the center of the Earth and the detector) 
are 1°.7 N, 0°.3 N, 61° . 7 N, and 73°.8 S, respectively. 
Following ICowan et al.l (|201l[ ). we henceforth refer to 
these 4 observations as Earthl:equinox, Earth5:solstice, 
Polar linorth, and Polar2:south. 

Figure [1] displays the 1-2/im portion of the observed 
reflection spectra of the Earth in Earthl. The broad 
absorption features at 1.08-1.18, 1.30-1.53, and 1.75- 
1.99/im are mostly due to H2O, and the narrower fea- 
ture around 1.27/im is due t o O2 plus oxygen collision 
complexes 02-02 and 02-N2 (jPalle et al.ll2009L and ref- 
erences therein). Abso rptions at 1.6/im and 2.0/im are 
signatures of CO2 (e.g. iRobinson et alll201lt ). 

We measure the equivalent widths of H2O (i^h2o) 
O2 ('^02) foi" each exposure. For wyl^q^ we focus on the 
spectral features centered at ~ 1.13/im and consider the 
absorption from 1.07/im to 1. 24/im. For wq^, we use the 
spectral features centered at ~ 1.27/im and consider the 
absorption from 1. 24/im to 1.283/im. In each case, the 
continuum line is assumed to connect the data points at 
both boundaries linearly. Additionally, we consider the 
variation of reflectivity at 1. 24/im as a measure of the 
continuum level. 

Figure [2] shows the diurnal fluctuations of wu^o and 
wq^ , as well as the variation of the reference continuum 
level (1.24 /im). Table [T] summarizes the average and 
the fractional variation amplitudes of continuum level, 
WYL2O and wq^. The variation amplitude defined by 
(maximum- minimum)/average is typically 5-20%. Fig- 
ure [2] clearly indicates that the diurnal patterns of wu^o 
and wq^ significantly differ from that of the continuum 
level which primarily traces the continental distribution 
on th e Earth (see e.g. Cowan et al. 2009 , 201li : lFujii et all 
l201lf ). Although not shown here in detail, we also con- 
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Earth 1: equinox 



Earths : solstice 



Polar 1: north 



Polar2: south 




Fig. 2. — Diurnal variations of continuum level (top), equivalent width of H2O absorption centered at 1.13/im (middle), and that of O2 
centered at 1.27/im (bottom). Symbols indicate observed data points and the connecting lines are drawn only to guide the eyes. Snapshots 
at corresponding times were generated at http: / / www.fourmilab.ch/cgi-bin/ uncgi/Earth are attached. 



TABLE 1 

Diurnal Variation of Continuum Level, Equivalent Width of H2O at 1.13/iM and that of O2 at 1.27/iM. 





Reflectivity at 1.24/im 
ave. (max-min) / ave 


H2O at 1.13/im 
ave. [A] (max-min)/ave 


O2 at 1.27/im 
ave. [A] (max-min)/ave 


Earthl:equinox (Mar. 2008) 
Earth5:solstice (Jun.2008) 
Polarlmorth (Mar. 2009) 
Polar2:south (Oct. 2009) 


0.111 32.1% 
0.084 30.6% 
0.084 29.5% 
0.084 18.2% 


457 16.5% 
499 18.2% 
392 16.1% 
372 9.4% 


39 16.8% 
42 18.7% 

46 7.2% 

47 5.5% 



firmed that absorption features of H2O at otiier wave- 
lengtiis (centered at ~ 1.4/im and ~ 1.85/im) exiiibit 
variation patterns matciiing tiiose at tiie 1.13/im band. 
For equatorial observations, the general trend that the 
depth is weaker in the first half of the day and becomes 
stronger in the second half is evident, and is shared by 
both H2O and O2. Referring to the snapshots shown 
in the bottom panels in Figure [21 the weaker absorption 
corresponds to the time when the Indonesia, a persis- 
tently cloudy region, do minates the field-of-view (see also 
iGomez-Leal et al.]l2012f ). This already demonstrates the 
strong relation between the absorption bands and cloud 
cover. In the next subsection, we investigate this con- 
nection in more detail. 

2.2. Correlation between diurnal variabilities and 
ocean/cloud parameters 

As mentioned in Section 12. If the variation pattern in 
Figure [2] is likely correlated with the extent and nature 
of cloud cover. In order to confirm the correlation be- 
tween the absorption features and clouds, we collect daily 
global maps of atmospheric parameters for the corre- 
sponding days from Terra/MODIS Atmosphere Level 3 
Product, which is available online^. Each global map 
is derived on a pixel-to-pixel basis from the Remote 
Sensing data obtained with MOderate Resolution Imag- 
ing Spectroradiometer (MODIS) o nboard the Ear th Ob- 
serving Satellites Terra and Aqua (jPorothv et al.l [2QQ6). 
Among the various parameters reported, we choose four 
which we suspect will influence the diurnal absorption 
variations strongly, including Atmospheric Water Vapor 

^ |http: / /ladsweb. nascom.nasa.gov/"] 



Mean^ Cloud Top Pressure Mean (Pctp), Cloud Frac- 
tion Mean (/cid), and Cloud Optical Thickness Combined 
Mean (tcm). We further add the ocean fraction as a fifth 
parameter. 

We compute the weighted average of each parameter for 
each EPOXI exposure. We adopt the geometric weight 
as Max{/io/^i, 0}, where /io/^i denote the cosine between 
the normal direction of the surface and the direction to- 
ward the Sun/observer. We adopt this factor because 
each surface pixel would contribute to the total reflec- 
tivity with that weight if the scatterings; were isotropic 
(Lambert's Law) (e.g. Les ter et al.lll979l ): while this is 
not strictly the case, it is a reasonable approximation for 
these purposes. 

Figure [3] displays Spearman's rank correlation coeffi- 
cients between absorption depths of II2O (left) and O2 
(right) measured by EPOXI observation and parameters 
taken from MODIS. Spearman's rank coefficient is cal- 
culated by converting each value Xi (z: index for exposures 
in one series of observation) to the rank Xi. If Xj is the 
n-th largest value of all possible values of x^, the rank 
is defined as Xj = n. Then we calculate the correlation 
coefficient of the ranks of two different parameters Xi and 
Vi as 

Y,{X^-X){Y,-Y) 
rs= ' — I (1) 

x=^x„ Y = Y,y^' (2) 

i i 

Since the chosen parameters do not necessarily follow a 
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Fig. 3. — Spearman's rank correlation coefficient between ab- 
sorption depth (left:H20, right :02) and Atmospheric Water Vapor 
(WV), Cloud Top Pressure (Pctp), Cloud Optical Thickness (rftp), 
Cloud Cover Fraction (/dd)? and Ocean Fraction (Ocean), respec- 
tively. Horizontal dotted lines show the significance level of 0.05 
for 13 samples (# of exposures per observation). 

gaussian distribution, we adopt Spearman's rank corre- 
lation coefficients instead of Pearson's correlation coeffi- 
cients. 

In most cases, the absorption depths are positively cor- 
related with the cloud top pressure (or equivalently, anti- 
correlated with the cloud altitude), while anti-correlated 
with both cloud optical thickness and the cloud cover 
fraction. These trends are plausibly understood as fol- 
lows: For cloud top pressure, the incident light is scat- 
tered at the altitude of the upper cloud layer without 
absorption by molecules in the atmosphere below the 
cloud. Since molecules are abundant in the lower atmo- 
sphere, the higher cloud altitude reduces the absorption 
depth more effectively. As for cloud optical thickness, 
thicker clouds scatter the incident light more completely 
and thus reduce absorptions below that layer. Addition- 
ally, the higher reflectivity of the thicker cloud increases 
the contribution of the region to the disk-integrated spec- 
tra. Similarly, the absorption depths become weaker as 
the cloud cover fraction increases. 

Although more weakly, the absorption depths appear 
to be slightly anti-correlated with the ocean fraction. 
The lower reflectivity of ocean reduces the absorption 
depths of the disk-integrated spectra^. 

We should emphasize here, however, that the above 

^ The light reflected from ocean in principle contains the sig- 
nature of liquid water (Palmer & Williams 1974), which is shifted 
toward longer wavelength and less prominent compared to the sig- 
nature of water vapor we consider in this paper. 



interpretation requires some caution because the five pa- 
rameters we adopted are not necessarily independent. In 
the equatorial region, for instance, there exists a positive 
correlation between cloud optical thickness and cloud top 
altitude. There is also generally positive correlation be- 
tween cloud cover and atmospheric water vapor, likely 
leading to the counter-intuitive anti-correlation between 
the atmospheric water vapor and the absorption depths. 

Figure [3] also shows that the correlation of the H2O 
band with clouds varies more significantly than that of 
O2, depending on the location of the sub-observer's lat- 
itude. In particular, the correlation of H2O for Polar 2 
observation even changes its sign compared to the other 
three cases. This is not the case for O2. We will discuss 
this behavior in detail in Section IH 



3. A SIMPLE OPAQUE CLOUD MODEL 

In the previous section, we discussed the influence of 
clouds on absorption features by considering the cor- 
relation between the absorption depths and individual 
cloud/surface parameter. In this section, we try to in- 
terpret the observed variation in absorption depths more 
quantitatively. For that purpose, we consider a simple 
opaque cloud model (see also e.g. iHarrison fc CoombesI 
Il988). which assumes that the cloud cover totally blocks 
the influence of the lower atmosphere and that the inci- 
dent light is completely scattered at the cloud top. 



3.1. Model Description 

The equivalent width of an atmospheric absorption 
band in the disk-integrated scattered light of the planet 
can be expressed as 



dfl = sin OdOd<p^ ^(^, 0) = /ioMi, 



(3) 

(4) 



where (^, (j)) represents the polar-coordinate of the sur- 
face point on the planet, /io = ^0(6^, ^) and /ii = (/)) 
are the directional cosines between the normal direction 
of the surface point and the direction toward the Sun and 
the observer, respectively. The wavelength-dependent 
continuum level (i.e., reflectivity) and optical depth are 
denoted by cxiO^cj)) and tx. The interval of integration 
over A is centered at the center of the absorption band, 
Aq. The integral over solid angle of planetary surface 
is performed over the illuminated and visible portion, 
Sisj. The above expression is based on the approxima- 
tions that 1) the lower boundary scatters the light ac- 
cording to the Lambert Law, i.e. the scattered intensity 
is independent of the emergent direction, and 2) scatter- 
ing by the atmosphere is negligible. 

We further neglect the wavelength-dependence of the 
continuum level, and assume that the wavelength depen- 
dence of cxiO^cj)) can be factored out, i.e., cx{0^(j)) = 
h{X)c{0^ (j)) with h{X) being a function of wavelength 
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alone. Then, equation (j3j) may be approximated as 



^model 



_ jg^^c{o,4>)w{N;e,4>)g{e,4>)dn 



(5) 



W{N; e,4>) = J d\ (1 - exp{-T^{N; 0, 4>)}) , (6) 

Tx{N; e, <P) = r n{z)ax {— + — ) dz 
Jzo \ Mo Ml / 



1 l ^ / ^{z)dz ■ J no{z)axdz 

Mo Ml/ / noiz)dz 

No 

1 1 



-'''A.O , 



A^eff=( — + — ) / n{z)dz, 
.Mo Mi/Ao 



(7) 
(8) 



where W{N;6^(j)) is the equivalent width at each sur- 
face patch, ax is the absorption coefficient, Nq is the 
total column density of molecules under the canonical 
atmospheric model, A^eff is the effective column density 
obtained by integrating n{z) along the optical path from 
the top of atmosphere down to the boundary at altitude 
z = zq^ and rA,o is the absorption depth of the canonical 
model. Equation ([7j) is only approximately valid because 
in reality the absorption coefficient ax depends on the 
T-P profile of atmosphere. Nevertheless, we employ this 
approximation for the sake of simplicity. 

In practice, we divide the planetary surface into 2° x 2° 
pixels, and treat the cloudless and cloudy portions in 
each pixel separately for convenience, as described below. 
Then, equation ([5]) is discretized as 



E 



{(1 



J eld/ ''^'^ nocld^nocld 



'^model 



E - ■^cld) Ccld + /eld 4ld} 9"^^" 



(9) 

where i is the index for the surface pixels, /dd is the cloud 
cover fraction, and the suffix nocid/cid indicates the value 
at the cloudless/cloudy portion. The position-dependent 
parameters in equation (|9j) are determined on the basis 
of daily global maps which are extracted from the same 
MODIS data product used in Section 12.21 For instance, 
the value for /dd is adopted from Cloud Fraction Mean. 
Input data for other parameters are described below. 

The effective column density A/eff of H2O at each patch 
is determined by Atmospheric Water Vapor Mean (w) 
and Cloud Top Pressure Mean (Pctp) as well as the ge- 
ometric factors /io and jui. For the cloudless portion, 
A/eff is identical to {1 / jno -\- 1 / jiii)w because we can safely 
set zq = 0. For cloudy portions of a pixel, however, we 
need to define the boundary height zq in equation ([7]) 
that is essentially the layer below which the the molecule 
in question does not influence the emergent spectra. In 
the present opaque cloud model, we assume that the syn- 
thetic absorption depth is completely unaffected by the 
atmosphere below the cloud top altitude Zctp, and thus 
set zq = Zctp. According to this opaque cloud assump- 
tion, N is determined once the vertical profile of water 
vapor and the altitude of the cloud top layer are given. In 
what follows, we simply assume that the vertical profile 
of water vapor is proportional to that of the US Standard 



Model, and determine the proportionality coefficient so 
that the total column density matches the value of w. 
The altitude of the cloud top layer is derived from Pctp 
assuming the US Standard Temperature-Pressure (T-P) 
profile^. The total column density A^eff for O2 is esti- 
mated in the same way except that we neglect the hor- 
izontal inhomogeneity, and determine it only from Pctp 
and the geometric factor. 
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Fig. 4. — Curve of growth for the equivalent width of the wa- 
ter vapor absorption band at 1.13/im-1.24/im, oxygen at 1.24/im- 
1.283/im, and carbon dioxide at 1588/im-1.622/im calculated by 
lbl2od with GEISA 2011 database under the no-cloud condition 
with surface albedo 0.3. 



The remaining task is to compute the equivalent width 
W as a function of the effective column density. For 
this purpose, we compute line-by-line optical depths rA,o 
under the US Standard Atmosphere with the GEISA 
2011 molecule spectroscopic database^^ and the pub- 



lic code lbl2od^^. Then we calculate /aq-aa 
with varying k,. Figure [4] shows the resultant curve of 
growth of the equivalent widths of II2O measured at 
1.06 - 1.24/im, O2 at 1.235 - 1.310/im as weh as CO2 at 
1.560 — 1.630/im used in Appendix [XI Based on these cal- 
culations, the normalized effective column density of each 
patch, A^gff/A/'usstandard, is translated into the equivalent 
width W by equating TV^ff/^uSstandard to in Figure H 
Finally, the continuum level c in equation Q is es- 
timated using the 2-stream approximation (e.g. iLioul 
1980). We consider a non-absorbing atmosphere with 
optical thickness and asymmetry factor p. Denoting 
the reflectivity at the surface by r^, the net reflectivity 
at the top of the atmosphere is: 



c(rcid) 



I -rn 



l + (V^/2)(l-r^)(l-/3)rcid 



, (10) 



where we adopt the optical thickness of the cloud ( Cloud 
Optical Thickness Combined Mean) for r^d- The typical 
value for /3 of Ear th's clouds at the relevant wavelengths 
is /3 = 0.85 (e.g. ILioul I1980D . The overah scaling fac- 
tor a is introduced here to empirically incorporate the 
anisotropic scattering due to the bulk cloud cover. We 

^ T-P profiles are not significantly different in different atmo- 
spheric models. 

http:/ /ether. ipsl.jussieu.fr/etherTypo/?id=1293 
http://www.libradtran.org/doku.php?id=lbl2od 
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determine the value of a by hand so as to match the 
disk-averaged continuum level to the observed data. 

According to the procedures described above, we com- 
pute VF* and at each pixel, and then obtain the pre- 
dicted absorption depth in the disk-averaged spectra of 
the Earth from equation (|9]). The model prediction is 
compared with the EPOXI data in Section [321 

3.2. Comparison of the opaque cloud model and the 
EPOXI Data 

Figure [5] compares the EPOXI data (symbols) to our 
model predictions (dotted lines). The left, middle and 
right panels correspond to the disk-integrated continuum 
level, H2O absorption, and O2 absorption, respectively. 
The anisotropic parameter for bulk cloud scattering a 
in equation (p!Q|) is determined so that the disk-averaged 
continuum level matches the observed value (Table [2j). 

TABLE 2 

Normalization Parameters and Comparison in Daily 
Average and Standard Deviation between Simulations and 
Observations. 





Earthl 


Earth5 


Polar 1 


Polar2 


a 


0.61 


0.64 


0.73 


0.79 


'"^H20, model/ '"^H20,obs 


0.85 


0.87 


0.88 


0.88 


'"^02, model /'^02,obs 


0.65 


0.66 


0.64 


0.63 


<^H20,model/f^H20,obs 


1.32 


1.73 


0.70 


0.92 


f^02,model/f^02,obs 


0.57 


0.77 


0.61 


0.88 



We first consider the comparison between the simula- 
tion and the observation in terms of the daily average, 
iD, and the standard deviation, a. Table [2] summarizes 
the ratio of the daily average of the simulated iDmodei, to 
the observed one, Wohs^ '^modei/'^obs, and that of stan- 
dard deviation amodei/^obs- The daily simulation average 
typically results in a 10-15% underestimation for H2O 
and 35% underestimation for O2. The underestimation is 
likely due to our "opaque cloud" assumption which tends 
to diminish the absorption depth; in reality cloud cover 
does not completely prevent the light from going through 
the lower atmosphere. Uncertainty can also come from 
the difficulty in determining the continuum level because 
the absorption bands we are considering are surrounded 
by other absorption features. 

The variation pattern, which is of our primary interest 
in this paper, is less sensitive to those uncertainties. The 
absorption depths normalized so that the daily average is 
and the standard deviation is unity are exhibited in the 
middle and right panels of Figure [5l In most cases, the 
overall variation patterns of both H2O and O2 are well re- 
produced by our opaque cloud model. In particular, the 
equatorial data (Earthl, Earth5) exhibit striking agree- 
ment with the model predictions. For polar observations, 
the agreement is somewhat degraded, especially for H2O 
variations in the Polar2 data. This may be partly as- 
cribed to the fact that there is a slight mismatch between 
the time of the EPOXI observation and that of the in- 
put atmospheric data of our simulation; the local time 
of MODIS observation is 10:30am for Terra and 12:10pm 
for Aqua^ while the disk-integrated spectra observed by 
EPOXI reflects the information of slices with different 
local times. In addition, the assumed vertical profile for 
H2O (the US standard atmosphere) is not expected to be 



as accurate for the equatorial/polar regions. While our 
model is admittedly rough in this regard, it is encour- 
aging that such a simple model reproduces the general 
trends and features observed by EPOXI fairly well. 

4. DIFFERENT VARIATION PATTERNS OF H2O AND O2 
ABSORPTION DEPTHS 

We now focus on the intriguing differences in behavior 
between H2O and O2. In the Earthl and Earth5 data, 
for instance, note the small deviation of the H2O pat- 
tern from the O2 pattern, including the sharper bump 
of W02 at t = 18 [hr] than il'h20- Polar 1 observations 
show the variation of k;h20 rising at 11 < t[hr] < 19, 
which is not present in the variation of W02. In addi- 
tion, only W02 has a local peak at t = 6. Significantly 
and reassuringly, these divergences are reproduced by our 
simulations, which allows us to isolate the origin of this 
difference. 

Let us discuss the possible causes of the different di- 
urnal pattern of II2O and O2 bands. Given that the ab- 
sorption depth in our model is directly linked to the total 
number of molecules above the cloud layers, the different 
behaviors of II2O and O2 exhibited in Figure [5] can only 
be ascribed to the differences in their column number 
density distribution in the planetary atmosphere. 

In the case of the Earth, the distribution of water vapor 
differs from that of O 2 in at least two ways. Firstly, the 
total column density of atmospheric water vapor is, un- 
like O2, a strong function of latitude (highly concentrated 
in the equatorial region) with additional local/temporal 
fluctuations. Secondly, the vertical profile of water vapor 
mixing ratio typically decreases as a function of altitude 
due to condensation, in contrast to the constant mix- 
ing ratio (volume fraction of the molecule) of O2 in the 
troposphere (Figure [6|). 

In order to evaluate these effects quantitatively, we run 
simulations of II2O band variation with different assump- 
tions for its atmospheric distribution. We consider two 
model distributions in addition to our fiducial model de- 
scribed in Section [STTl Model A adopts the water vapor 
vertical profile of the US Standard Model everywhere, 
and thus the mixing ratio decreases with altitude (red 
solid line in Figure [6]). Model B instead assumes a con- 
stant mixing ratio of H2O in the troposphere as is the 
case for O2 (blue dashed line in Figure [6]). Both Models 
A and B neglect the inhomogeneity of water vapor over 
the surface of the planet. 

These three simulations and observations in the Po- 
lar l:north case are plotted in Figure [71 The observed 
variation of II2O agrees well with fiducial model, while 
the other two models that neglect the horizontal or sur- 
face inhomogeneity and/or peculiar vertical profiles fail 
to match the observed data. We also note that the Model 
B (blue line in Figure [7j) prediction of wn^o is very close 
to that of the simulation for O2 in the Polarl case (solid 
line in the third panel of O2 in Figure [5]). This indeed 
implies that that the differences in variation patterns ba- 
sically reflect their differing spatial distributions in the 
atmosphere. 

Both of the above-mentioned properties of water vapor 
distribution are consequences of the coexistence of mul- 
tiple phases of water on Earth. The horizontal inhomo- 
geneity of water vapor originates from the temperature 
gradient over the surface (since the vaporization rate is 
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Fig. 5. — Comparison between the EPOXI data (symbols) and simulations (dashed lines) in terms of continuum level measured at 1.24/im 
(left), water vapor absorption (middle), and oxygen absorption (right). 
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Fig. 6. — Vertical profiles of H2O (red solid) and O2 (blue 
dashed) taken from US standard model. 

basically determined by the surface temperature) and/or 
the inhomogeneous distribution of liquid/solid reservoirs 
on the surface. The water vapor distribution is not mixed 
well since the mean residence time (MRT) of atmospheric 
H2O (~10 days) is much shorter than the time scale of 
the vertical and global mixing times of the atmosphere 
(both - 1 year). In contrast, the MRT of O2 is -4000 
years, thus insuring that O2 is well mixed and only very 
slowly changing, if at all. 

The short MRT of water vapor implies that there is a 
non- negligible flux of H2O into/out of atmosphere and 
the existence of a substantial water reservoir (compared 
to the abundance in the atmosphere), such as the oceans 
on the Earth. Figure [7] also indicates that the difference 
in the vertical profile alone may lead to differences in the 
variation of absorption bands, but this is also related to 
the condensation of water in the atmosphere. In either 



Fig. 7. — Comparison in H2O band variation between obser- 
vation and simulations with different assumption for water vapor 
distribution. Model A adopts the water vapor vertical profile of US 
standard model but neglect the horizontal inhomogeneity. Model 
B assumes constant mixing ratio in the troposphere as is the case 
for O2. 

case, the deviation of absorption features of H2O from 
the well-mixed gas like O2 may serve as an indicator of 
phases transitions happening in the atmosphere or on the 
surface of a planet. This signature of the Earth's "wa- 
ter cycle" could be an observable indicator of a similar 
climate on a terrestrial exoplanet. 

5. SUMMARY AND DISCUSSION 

In this paper we examined the diurnal variability of 
the absorption depths of the two important molecules, 
H2O and O2, in the disk-integrated spectra of the Earth. 
We analyzed the H2O band centered at ~ 1.13/im and 
O2 band centered at ~ 1.27/im observed by EPOXI and 
found that these absorption bands show diurnal vari- 
ations correlated with uneven cloud cover. A simple 
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opaque cloud model, which assumes that the cloud com- 
pletely blocks the atmospheric signatures below the cloud 
top layer, is able to reproduce the basic variation patterns 
of the H2O and O2 bands using pixel- level cloud data ob- 
tained with Earth Observing Satellite. Thus we conclude 
that the non-uniform cloud cover distribution dominates 
the observed diurnal variations of those molecular ab- 
sorption depths. 

However, we also found that the diurnal variability pat- 
terns of H2O and O2 bands are not identical, and the dif- 
ferences originate from the inhomogeneous distribution 
of water vapor in the atmosphere. The variability pat- 
tern of O2 is basically explained by the the cloud cover 
distribution because it is well mixed in the troposphere 
of the Earth. In contrast, the variability pattern of H2O 
is well reproduced only with additional information on 
the vertical profile and spatial inhomogeneity of atmo- 
spheric water vapor as a model input. The nature of 
the water vapor distribution in the atmosphere is linked 
to the fact that H2O circulates in the planetary surface 
layer changing its phase among water vapor, liquid wa- 
ter (ocean/lake/pond/river) and/or water ice. There- 
fore, different behavior in the variability patterns of O2 
and H2O absorptions may carry information on the in- 
homogeneous phase distribution of H2O in the surface 
layer. 

Our study of the Earth demonstrates the possible role 
of the variability of absorption bands in characterizing 
the surface environments of Earth- analogs. If future in- 
struments eventually succeed in detecting scattered light 
from such exoplanets, their time-resolved measurements 
of the absorption bands will reveal the uneven cloud 
cover and/or inhomogeneous spatial distribution of at- 
mospheric constituents. Our current results suggest that 
the comparison in variation pattern between presumably 
well-mixed gases (such as O2) and II2O may indicate 
frequent phase transitions of water in the surface layer 
and thus serve as a probe of habitability, in a comple- 
mentary fashion to the direct detection of liquid wa- 
ter (e.g. IWilliams GaidosI [2008[ lOaklev C^[2009[ 



iRobinson et al .1120101 : IZugger et al.ll201Q[) and the search 
for other potential bio si^natures (e.g. iKaltenegger et al.l 
I2OIQI : iKawahara et al.l [20 12 ) . 

While we focused on absorption bands in the NIR 
range where the EPOXI data are available, it is natural 
to expect that other bands also exhibit similar variation 
patterns at wavelengths dominated by scattering of the 
primary star's light, rather than planetary thermal emis- 
sion. In the visible range, there are O2 bands at 0.69/xm 
(equivalent width ~13.5A) and at 0.76/im (equivalent 
width -47.8A) fe.g. lPalle et~an[2QQ9. ) as weh as several 
H2O bands. 

Inhomogeneity of an atmospheric constituent becomes 
appreciable when its MRT is short compared to the time- 
scale of the atmosphere mixing. This is equivalent to 
there being a substantial flux of the molecular species 
in question between the atmosphere and some surface 
reservoir, i.e., substantial compared to the total amount 
residing in the atmosphere. Besides vaporization, photo- 
chemical production (e.g. O3) or biotic production (e.g. 
N2O, CII4) could lead to such a situation. If sufficiently 
accurate data become available in the future, we may be 
able to interpret implied inhomogeneous distributions in 
terms of those or other specific mechanisms. 
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APPENDIX 
VARIATION OF CO2 

We also examine the variability of CO2, another absorption feature imprinted in the observed spectra. Absorption 
bands of CO2 exist around 1.6/im and 2/im (Figure [T]). These CO2 bands override the wings of H2O absorption, 
making it harder to determine the continuum level. We pick up the absorption band at 1.59-1.62/im, which is likely to 
be least affected by H2O absorption. The peak-to-throat variation amplitude of four observations are 5-20%^^ (Table 
[3]), on the same level of O2 and H2O. 

Figure [8] shows the diurnal variation patterns of CO2 band at 1.59-1.62/im extracted from EPOXI data as well as 
those of our simulation following the same procedure as H2O/O2 and assuming the uniform distribution (the assumed 
growth curve of this band is displayed in Figure H]). The normalization factors are summarized in Table HI We see the 
clear similarity in variation patterns to O2 rather than H2O. For instance, the clearer bumps at t = 18 of Earthl and 
Earths observations and the gradient in Polar2 observation. This result is consistent with our discussion in Section SI 
since the MRT of CO2 is 3-5 years and slightly longer than the tropospheric mixing timescale. The behavior of Polarl 
is less conclusive and the observation does not match the simulation so well . Although we do not identify the cause, 
it might be related to some level of inhomogeneity of CO2 due to the intermediate MRT. 

TABLE 3 

Diurnal Variation of Continuum Level, Equivalent Width of CO2 at 1.6/iM. 





CO 


2 at 1.6/im 




ave.[A] 


(max-min) /ave 


Earthl :equinox (Mar. 2008) 


12.1 


80.1% 


Earth5:solstice (Jun.2008) 


17.8 


17.4% 


Polarlmorth (Mar. 2009) 


17.4 


11.3% 


Polar2:south (Oct. 2009) 


17.9 


5.4% 



Carbon Dioxide (L59-L62)im) 




6 12 18 

Time [hr] 



Fig. 8. — Diurnal Variations of CO2 absorption bands around 1.59-1.62 /im. EPOXI data (symbols) and simulations (dashed lines) are 
compared. 



The large peak-to-throat amplitude of Earthl :equinox is at least partly due to the failure in continuum determination. In 
Earthl :equinox only, the observed wavelength grid changes for each exposure and the sharpness of the very narrow feature at 1.58/im, 
which is critical to determine the continuum level, also changes. This results in the uncertainty in estimation of equivalent width. 
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TABLE 4 

Normalization Parameters and Comparison in Daily Average and Standard Deviation between Simulations and 

Observations. 





Earthl 


Earths 


Polar 1 


Polar2 


^C02,model/'"^C02,obs 


1.14 


0.85 


0.95 


0.90 


C02,model/<^ C02,obs 


0.23 


1.23 


0.74 


1.75 



